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尼 米 的 引入 扩展 了 传统 厄 米 量子 系统 中 的 概念 并 诱导 出 许多 新 奇 的 物理 现象 ， 比 如 非 厄 米 系统 所 

独 有 的 非 厄 米 趋 肤 效 应 , 这 使 得 对 非 厄 米 量 子 模 型 的 模拟 成 为 大 家 关注 的 热点 . 相 比 于 量子 平台 , 经 典 

系统 具有 成 本 低廉 、 技 术 成 熟 、 室 温 条 件 等 优势 , 而 其 中 的 经 典 电路 系统 则 更 加 灵活 , 原则 上 可 以 模 卸 

任意 维度 、 任 意 格 点 间 跃 迁 、 任 意 边界 条 件 下 的 量子 紧 束缚 模型 , 已 经 成 为 模拟 量子 物 态 的 有 力 平台 . 
pg 


本 文 利用 经 典 电 路 通过 SPICE 成 功 模拟 了 一 个 重要 的 非 厄 米 量子 模型 非 互 易 Aubry-André 模 
型 — 的 稳 态 性 质 , 此 模型 同时 具有 非 互 易 的 格 点 跃迁 和 准 周期 的 格 点 在 位 势 . 以 此 为 例 , 详细 介绍 


如 何 建立 经 典 电路 的 拉 普 拉 泣 形式 与 量子 紧 束 缚 模型 哈密 顿 矩阵 在 不 同 边界 条 件 下 的 映射 ， 尤其 是 如 
何 利用 电流 型 负 阻抗 变换 器 构建 模型 的 非 互 易 性 , 然后 , 根据 电路 的 格林 函数 , 通过 AC 电流 驱动 并 六 
量 电压 响应 的 方式 , 用 SPICE 模拟 了 周期 边界 条 件 下 的 复 能 谱 和 相应 的 能 谱 绰 绕 数 ， 以 及 开 边 界 条 从 

的 趋 肤 与 局 域 模式 的 竞争 . 其 中 , 为 了 使 电路 的 响应 不 发 散 , 本 文 还 解析 地 给 出 辅助 元 件 的 设置 原则 . 
结果 显示 , SPICE 模拟 与 理论 计算 很 好 地 吻合 , 为 进一步 的 实验 实现 提供 了 详细 的 指导 ， 由 于 本 文 
路 设计 与 测量 方案 的 普 适 性 , 原则 上 可 以 直接 应 用 于 其 他 非 厄 米 量子 模型 的 电路 模拟 . 
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The introduction of non-Hermiticity into traditional Hermitian quantum systems generalizes their 
basic notions and brings about many novel phenomena, e.g., the non-Hermitian skin effect that is 
exclusive to non-Hermitian systems, attracting enormous attention from almost all branches of 
physics. Contrary to the quantum platforms, classical systems have the advantages of low cost and 
mature techniques under room temperature. Among them, the classical electrical circuits are more 
flexible on simulating quantum tight-binding models in principle with any range of hopping under 
any boundary conditions in any dimension, and have become a powerful platform for the simulation 
of quantum matters. In this paper, by constructing an electrical circuit, we simulate by SPICE the 
static properties of a prototypical non-Hermitian model 一 the nonreciprocal Aubry-André (AA) 
model that has the nonreciprocal hopping and on-site quasiperiodic potentials. 

The paper is organized as follows: Following the introduction, in Sec. II we review in detail the 
Laplacian formalism of electrical circuits and the mapping to the quantum tight-binding model. 
Then, in Sec. III, an electrical circuit is proposed with resistors, capacitors, inductors, and the 
negative impedance converters with current inversion (INICs), establishing a mapping between 
the circuit’s Laplacian and the non-reciprocal AA model’s Hamiltonian under periodic boundary 
conditions (PBCs) or open boundary conditions (OBCs). Especially, the nonreciprocity, the key of 
this model, is realized by INICs. In Sec IV, based on the mapping, for the proposed circuit under 
PBCs, we reconstruct the circuit’s Laplacian via SPICE by measuring voltage responses of an AC 


current input at each node. The complex spectrum and its winding number v can be calculated 


by the measured Laplacian, which are consistent with the theoretical prediction, showing v = +1 


for non-Hermitian topological regimes with complex eigenenergies and extended eigenstates, and 
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v = 0 for topologically trivial regimes with real eigenenergies and localized eigenstates. In Sec V, 
for the circuit under OBCs, a similar method is used for measuring the node distribution of the 
voltage response, which simulates the competition of non-Hermitian skin effects and the Anderson 
localization, depending on the strength of quasiperiodic potentials; the phase transition points also 
appear in the inverse participation ratios of the voltage responses. 

During the design process, the parameters of auxiliary resistors and capacitors are evaluated for 
obtaining stable responses, because the complex eigenfrequecies of the circuits are inevitable under 
PBCs. Our detailed scheme can directly instruct further potential experiments, and the designing 
method of the electrical circuit is universal and can in principle be applied to the simulation for 


other quantum tight-binding models. 
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本 文 的 主要 目的 是 通过 对 非 互 易 AA 量子 模型 的 经 典 
子 紧 束缚 模型 , 便于 感 兴趣 的 


路 模拟 , 详细 介绍 如 何 用 电路 的 拉 普 拉 半 形式 模拟 
量 读者 利用 类 似 方法 模拟 其 他 量子 模型 , 以 及 为 实验 实现 提供 详细 指引 . 剩 下 的 内 
容 安 排 如 下 : 第 二 节 详 细 介绍 如 何 构建 经 典 电路 的 拉 普 拉 半 形式 与 量子 紧 束 缚 模型 的 映射 , 第 三 广 具 体 给 出 实现 
不 同 边界 条 件 下 非 互 易 AA 模型 的 电路 设计 方案 , 第 四 节 和 第 五 节 利 用 SPICE 分 别 模拟 非 互 易 AA 模型 在 周期 
边 下 的 能 谱 和 能 量 缠绕 数 以 及 在 开 边 界 条 件 下 趋 肤 与 局 域 模式 的 竞争 , 最 后 一 节 进 行 总 结 . 


IT， 经 典 电 路 的 拉 普 拉 辛 形式 与 紧 束缚 模型 的 对 应 


别 对 应 电路 的 连接 点 和 元 件 (9) 


任意 经 典 电路 组 成 的 网 络 都 可 以 用 一 个 图 (graph) 来 表示 ， 其 节点 和 边 分 


ba, bo). 比如 由 电阻 、 电 感 和 电容 (简称 RLC) 等 被 动 元 件 组 成 的 电路 , 元 件 各 自 的 物性 方程 为 
dI dV 
V=RI, V=L CF =L (3) 


ZEW VA 分 别 表示 元 件 两 端的 电压 差 和 通过 元 件 的 电流 ，( 忌 , L,C) 


分 别 为 元 件 的 电阻 、 电 感 和 电容 . 根据 


HEIKER (Kirchhoff) 电流 定律 , 利用 以 上 物性 方程 , 可 以 得 到 图 中 每 个 节点 关于 时 间 # 的 微分 方程 : 
d,, ,ad .时 i, 
d 
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其 中 , T(t) 和 Vat) Pa Re HER A n 的 外 界 输入 电流 和 对 地 电压 , 这 里 用 (Ram, LnmsCnm) 分 别 表示 从 节 


Fin 到 节点 m (下 标 9 代表 接地 ) 的 等 效 电阻 、 等 效 电感 和 等 效 电容 ,以 便 描 述 更 一 般 的 具有 非 互 易 特 性 的 元 
件 ， 通常 的 被 动 RLC 元 件 为 互 易 的 ， Bp Rym = Rmn = R, Lam = Linn = L, Chm = Umn = C. 上 述 方程 可 以 写 
成 更 紧凑 的 矩阵 形式 : 
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其 中 , L(t) 和 VE) 分 别 表 示 市 点 输入 电流 和 对 地 A (R,L,C) 分 别 为 等 效 电阻 、 等 效 电感 和 等 效 电 
容 构成 的 系数 和 矩阵, 其 矩阵 元 为 


Cun = Cng RE `> Cranes 


+ LV(t), (5) 
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m#n m#n mm 天 nm 
Cam = -Cnm (nÆM), Rem = -Ri (n¥#m), Lam = -Lm (W#M). (6) 
对 于 微分 方程 B= 考虑 具有 固定 频率 w 的 AC 电流 源 Z(t) 及 其 电压 响应 V(t), 其 形式 为 
T(t) =Te*, V(t) = Ve. (7) 
将 它们 代入 方程 E), 得 到 不 含 时 的 矩阵 方程 : 
= (iw +R+L/iw)V = J (w)V. (8) 
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依赖 于 驱动 频率 w. 如 果 没有 外 界 电流 输入 , 即 T(w)V 


“(w)Z = G(w)T, 


= 0, 则 det J (we) = 0 决定 了 电路 的 本 征 


(9) 
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其 中 , 9(w) = J! (w) 被 称 为 电路 的 格林 函数 (circuit Green’s function), 具有 阻抗 的 量 纲 . 

实际 上 , 对 于 具有 形式 (A) 的 稳 态 解 , 任意 电路 网 络 都 可 以 表示 为 拉 普 拉 注 的 形式 ， 比 如 包含 放大 器 的 主动 
路 fd, bo, 53 和 具有 非 线性 元 件 的 非 线性 电路 j6q) b1) 等 
为 了 利用 经 典 电路 系统 模拟 量子 紧 束缚 模型 , 可 以 将 电路 拉 普 拉 辛 矩阵 .7(w) 直接 与 紧 束 缚 模型 在 实 空间 
的 哈密 顿 量 矩 阵 H 相对 应 5. HAERE 了 (w) 的 本 征 方程 


T (wb? (w) = jn(w)W(W), TWD (o) = Fw) (w), (10) 


HUERTA ST EA EA EASE REPS ER. 这 里 ，j(w) AEE SS n 个 本 征 值 , 构成 本 征 导 纳 谱 ， 
可 以 完整 模拟 H 的 能 谱 . 特别 地 , 由 det I (we) = 0 可 知 , 本 征 频率 w 使 至 少 一 个 本 征 导 纳 为 零 , 即 jn (we) = 0. 
于 了 (w) 一 般 为 非 厄 米 和 矩阵 , 即 了 i(w) A T(w), 相应 的 本 征 矢 通常 包含 本 征 右 矢 和 本 征 左 矢 wy (w), 利用 电 
路 的 交流 分 析 , 可 以 得 到 了 (w) 的 右 本 征 模式 LL (w), 从 而 模拟 H 的 右 本 征 态 ; 左 本 征 态 可 以 用 Tt (w) 模拟 . 
通过 元 件 以 及 驱动 频率 的 设计 和 调节 , 了 (w) 具有 高 度 可 控 性 , 原则 上 可 以 模拟 任意 维度 、 任意 边 界 条 件 、 非 线性 
和 非 厄 米 等 非常 广泛 的 量子 模型 的 稳 态 性 质 . 

接 下 来 , 应 用 以 上 电路 的 拉 普 拉 辛 形式 构建 7 了 (w) 与 非 互 易 AA 模型 哈 窗 顿 量 和 矩阵 H 的 对 应 关系 , 从 而 对 
其 稳 态 性 质 进 行 模拟 , 包括 周期 边界 条 件 下 的 能 谱 和 缠绕 数 以 及 开 边 界 条 件 下 趋 肤 和 局 域 模式 的 竞争 . 


THT. 非 互 易 AA 模型 的 电路 拉 普 拉 辛 


在 交流 驱动 下 , 被 动 元 件 往往 呈现 出 互 易 性 , 这 是 由 最 基本 的 基 尔 霍 夫 电流 定律 决定 的 , 例如 电容 与 电感 的 
导 纳 Jo(w) = iwC 和 Jp (w) = 1/iwL 均 不 依赖 于 正 向 或 反 向 测量 即 可 表征 。 而 根据 公式 S) 的 描述 ， 要 实现 拉 
EMER SER AE Imaw) A Inm(w), 则 需要 电路 中 某 一 元 件 的 导 纳 值 依赖 于 测量 的 方式 ， 这 通常 需要 引 
入 主动 元 件 , 比如 电流 型 负 阻抗 变换 器 (negative impedance converter with current inversion, 简称 INIC) b3). 


如 图 lil) 所 示 的 INIC 由 放大 器 和 若干 线性 元 件 构成 , 根据 基 尔 霍 夫 电流 定律 , 容易 得 出 两 端的 输入 电流 分 别 为 
WOZ. 
n=“), =iwCi(V - Vi). (11) 
这 表明 一 般 情况 下 INIC 两 端 不 同方 向 的 导 纳 不 相等 : 
J1(w) = Bele Jy (w) = iwCr. (12) 


Z+ ” 
为 方便 起 见 , 理论 上 选取 INIC 中 的 阻抗 满足 24 = Z_, 使 得 元 件 两 端 流 问 放 大 器 的 电流 大 小 相同 方向 相反 , 即 
五 三 五 , 从 而 得 到 两 个 方向 符号 相反 的 导 纳 : J1(w) = —J,(w). 

另外 , 这 里 会 用 到 具有 负 值 的 元 件 (比如 负电 阻 等 ) 其 两 种 实现 形式 如 图 lil) 所 示 kgl, 它们 通过 放大 器 分 
别 实现 了 对 地 单 端 口 和 上 自由 两 端口 的 等 效 负 阻抗 ( 导 纳 )， 根据 基 尔 霍 夫 定 律 , 可 得 对 地 单 端 口 电路 [图 tke) 左 
图 ] 的 输入 阻抗 为 


Vy Vy 
mI Tee =“ 
类 似 地 , 自由 端口 电路 [E fe) 右 图 ] 两 端的 输入 阻抗 分 别 为 
Vi- Vj Vi- Vj V-V _ Vj- V; 
mo- mu ~ aaa morae S 


BI, Zji = Zij = —Z. 
利用 以 上 关键 元 件 , 在 由 电感 Lo 和 电容 Co 组 成 的 左手 传输 线 离 散 模型 (lumped-element circuit model for 
a left-handed transmission line) 的 基础 上 , 设计 了 如 图 lia) 所 示 的 非 互 易 AA 模型 的 电路 模拟 示意 图 , 由 


图 1. (a) 上 图 : 非 互 易 AA 模型 的 电路 模拟 示意 图 , 包含 N 个 有 效 电压 节点 Va (n = 1,… , N), 节点 间 元 件 Co 和 INIC(,) 模拟 
格 点 间 的 耦合 , 其 中 INCo 用 于 实现 关键 的 非 互 易 耦合 , 其 定义 见 (b); 接地 元 件 ne 模拟 格 点 的 在 位 势 ; X 相关 模 
块 和 开关 控制 边界 条 件 的 模拟 . 下 图 : X 相关 模块 的 定义 . (b) INIC 元 件 的 内 部 电路 图 , 由 理想 放大 器 、 阻 抗 Z4 和 目标 元 件 Cr 


(没有 Xo) 构成 , 可 以 实现 Vi, 两 端 不 同方 向 的 导 纳 不 同 ; INTC, 仅 需 将 INIC 中 的 目标 元 件 Cr 再 并 联 一 个 Xe 即 可 . (c) 负 阻 抗 
模块 pA. 左右 图 分 别 实现 对 地 单 端口 和 自由 两 端口 间 的 等 效 负 阻 抗 一 2, 其 中 理想 放大 器 上 的 标记 表示 输出 电压 与 输入 电压 的 关 
系 . 各 元 件 的 具体 功能 描述 详 见 正 文 . 


Fig. 1. (a) Upper panel: Sketch of an electrical circuit simulating the nonreciprocal AA model. It includes N voltage nodes 


+ 


Vn (n = 1,---,N) with elements Co and INIC(p) simulating the intersite couplings, where INIC(,) defined in (b) is the key 
element to realize the nonrecprocity, and the grounded elements (Lo, Ro, Co,r,n) simulating the on-site potentials; X modules 
and the switches control the simulation of boundary conditions. Lower panel: Definitions of X modules. (b) The internal 


circuit of the INIC, constructed by the ideal operational amplifier (opamp), impedance elements Z+, and the targeted element 


Cr (without X»), which can realize unequal effective input inductances from the two different ports V;,,; INIC, is defined by 
adding an extra X» module in parallel with Cr in INIC. (c) Modules of negative impedance pg. The internal circuits of the 
grounded one-port and the floated two-port negative impedances —Z for the left and right panels, respectively, where the 
labels of the ideal opamps represent the relation of the output voltage to the input voltages. See relevant texts for the detailed 


description of each element. 


RLC 和 INIC 等 元 件 组 成 , 包含 有 N MARCHE TI A, 对 应 模拟 模型 的 NN AEA. 根据 方程 全， 容易 写 出 此 电 
路 的 拉 普 拉 演 矩阵 : 


一 C1 Co 一 C7 0 XR 
Co 二 CT —-Cy Co 一 Cr 
. ， et 1 
J (w) = 一 化 ) 区 eu CT) T iwLp T ra 
Co 二 CT —Cn-1 Co—Cr 
Xr 0 Co 十 CT 一 CN 
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其 中 ， 
XR/L = CoCr- (w? La)! + (w? Le) T (iwRa)* + (iwRy)*, (16) 
HAE Cnr ER Loos 电阻 Rae 和 Ro ER PEX, E 为 单位 矩阵 ; 边界 条 件 由 图 La) 中 X 相关 模块 与 开关 


控制 : 当 两 端的 开关 同时 接 入 端口 0 并 且 X 相关 模块 (Xow 和 Xur) 均 开路 时 , 系统 对 应 于 开 边 
的 开关 同时 接 入 端口 

将 此 电路 的 拉 普 拉 辛 矩阵 (ia) Saks AA 模型 的 哈密 顿 量 甜 阵 Q) 相 比较 ， 可 以 建立 除 
矩阵  ) 外 两 者 间 的 映射 A eo H, 并 利用 对 应 的 无 量 纲 化 参数 得 到 如 下 等 式 关系 . 

对 于 电路 的 主体 部 分 (不 包含 边界 ), 由 


Co 十 CT aa -Cn _ 2A cos(2r 8n) 
Gea °° * Gea er 
可 得 
Cr Ch 2A cos(2r Gn) 
FS eh = 
Co ae Co cosha ’ 


这 里 将 Co 作为 参考 电容 . 由 此 可 以 理解 电路 主体 各 元 件 的 作用 : (Co, Cr, Lo) 构成 紧 束缚 模型 的 


2 时 , 下 文中 将 会 看 到 , 通过 X 相关 模块 参数 的 调节 可 以 模拟 具有 磁 通 的 周期 边界 条 件 . 


界 条 件 ; 当 两 端 


正比 于 单位 


D ( 


(17) 


(18) 


整体 参考 势 (了 


部 分 ) , Co 还 承担 格 氮 间 互 易 耦 合 的 作用 , 非 互 易 耦 合 和 变化 的 在 位 势 由 INIC 中 的 


别 实现 ; 随后 可 以 看 到 , 电阻 Ro 和 电容 C; 的 引入 是 为 了 使 
体 平移 . 需 注意 的 是 , Cn 随 着 市 点 n 的 变化 会 被 要 求 为 负数 ， 


对 于 边界 部 分 , 同样 利用 对 应 的 无 量 纲 化 参数 关系 
XR 二 Xr ee? t etet? Xr- Xr ee? — ete? 
2Co epee ’ 2C; eeo ’? 
可 得 
La we Ry Woy 
= 人 
Lo  w?(1—cos®)’ ° i Ro wsin® ii 


BA Cr 和 接地 电容 Cn 分 
路 的 响应 不 发 散 , 它们 仅 使 导 纳 谱 在 复 平面 内 作 整 
等 效 负电 容 可 以 使 用 图 也) 的 方案 实 


现 . 


(19) 


(20) 


这 里 将 Lo 和 Ro 分 别 作 为 参考 电感 和 参考 电阻 , wo = 1/V LoCo 作为 参考 频率 , 并 定义 无 量 纲 量 y = 
当 D 为 2r 的 整数 倍 时 , Lao 和 Ras 均 发 散 ， 代 表 开 路 ， 对 应 于 无 磁 通 的 周期 边界 条 件 . 随 
会 为 负数 , 同样 可 以 使 用 图 Mo 方案 实现 . 


TV. 周期 边界 条 件 下 能 谱 和 缠绕 数 的 模拟 


众所周知 , 即使 在 周期 边界 条 件 下 , 无 序 系统 也 不 再 具有 平移 不 变性 ， 
空间 的 方法 计算 系统 的 缠绕 数 , 而 通常 的 办 法 是 在 链 环 中 心 加 入 磁 通 量 为 更 的 磁场 , 此 时 系统 变 


因而 无 法 通过 将 哈密 顿 量 


Fg i bag 
P 的 变 化 ， Rap 


变换 到 动量 
为 ® 的 周期 函 


数 (周期 为 2r), 从 而 进行 计算 . 对 于 非 厄 米 系统 , 由 于 能 量 一 般 为 复数 , 可 以 定义 复 能 量 在 复 平 


四 的 缠绕 数 来 刻 


画 非 厄 米 系统 的 拓扑 相 7 ba: 


27 27 
-| Os In det H(® ) Jð = f Os o(d) 6 
0 0 


271 2T 


其 中 , 0(B) 是 det H(®) 的 幅 角 . 对 
其 边界 条 件 下 , v = 0 表示 拓扑 平庸 的 局 域 相 , v = 士 1 表示 两 种 拓扑 非 : 
为 了 通过 电路 模拟 并 测量 非 互 易 AA 模型 的 缠绕 数 v, 利用 两 者 的 对 应 关系 可 以 将 定义 (2) 


a? 


青 的 扩展 相 . 


~ 
=> 


FIRE AA 模型 , 由 文献 PA 可 知 , 不 同 的 缠绕 数 表示 不 同 的 拓扑 相 : 在 周 


(21) 


的 H 用 拉 普 


PFE RE (5) 的 A/Co 替换 (这 里 除 以 Co 是 为 了 保证 In 等 函数 的 作用 对 象 是 无 量 纲 的 ， 
影响 v 的 结果 ). 因此 , 只 要 能 从 实验 上 测量 出 不 同 于 下 的 A(®) 矩阵 , 即 可 计 


整体 的 倍数 并 不 会 
出 相应 的 缠绕 数 ， 


的 909 960 的 的 


Imlwe]/wo 
Im[an] 


-1 a = - 
Re[we]/wo 1.35 -6 Refan] 6 
(Dots (d) 1 * 
a > 
© a Lg 
ee Fe 6 < > 
= . s ° a =) Ree tee cee eee eee ee cee ree ce 
wile e ER = g , 
一 eee yA @ [ R a 
96600 Me AT 
0 一- 一 -人 一 一 一 一 -1 : - 4 ay 
0.32 Relwj/uo 052 -2 Refan] 2 0 b/r 2 


图 2. (a) 非 互 易 AA 模型 的 参考 相 图 , 相 边 界 (红色 实 线 ) 由 文献 j 附录 中 的 解析 表达 式 (B6) 计算 缠绕 数 v 求 得 . 虚线 用 了 
图 B 的 计算 . (b,c) 分 别 对 应 (a) 中 B、C 两 点 的 电路 本 征 频 率 谱 ( 取 oy = 0 时 ) ( 左 图 ) 以 及 矩阵 A/Co 的 本 征 谱 ( 取 & = 0 FY) 

图 ). 红色 实 点 和 蓝 色 实 方块 分 别 代表 周期 和 开 边 界 条 件 下 的 理论 结果 . 左 图 中 的 空心 菱形 标记 驱动 频率 w 的 位 置 , 右 图 中 的 
心 菱形 表示 周期 边界 条 件 下 SPICE 的 模拟 结果 ，(d) 根据 SPICE 模拟 得 到 的 A(©) 计算 得 出 0(®) = In det [4(@)/Co] BE 2 
变化 . 左 向 三 角 、 右 向 三 角 和 菱形 分 别 对 应 (a) 中 A、B 和 C 三 点 , 虚线 为 相应 的 理论 值 . 以 上 所 有 计算 取 有 限 尺 寸 N = 21, 电 
路 元 件 的 具体 设置 参见 相关 正文 . 
Fig. 2. (a) Referenced phase diagram of the non-reciprocal AA model, where phase boundaries (solid red lines) are obtained 
by calculating the winding number v with Eq. (B6) in the Appendix of Ref. Eh. The dashed lines are for the plots in 
Fig. 3. (b,c) The eigenfrequencies of the circuit at y = 0 (left panels) and the eigenvalues of A/Co at ® = 0 (right panels). 
Solid red dots and solid blue squares represents the theoretical results under PBCs and OBCs, respectively. The hollow 


EZA 


diamonds in left panels label the driving frequency, and those in right panels are the simulated results under PBCs via SPICE. 
(d) 0(®) = Indet [A(®)/Co] versus ® based on the simulated A(®) via SPICE, where left-pointing triangles, right-pointing 
triangles, and solid diamonds represent points A, B, and C in (a), respectively. The dashed lines are the corresponding 


theoretical results. All figures are calculated in a finite size N = 21. See relevant texts for the specific setting of the circuit’s 


elements. 


本 文 利用 电路 的 格林 函数 形式 @) 进行 SPICE 模拟 . 对 于 具有 周期 边界 (将 图 目 所 示 电 路 两 端的 开关 均 置 
到 p 端口 ) 的 电路 , 仅 在 第 n 节点 楼 入 频率 为 w 的 AC 电流 源 , 测量 所 有 N 个 节点 的 电压 响应 , 并 除 以 输入 
流 强度 , 即 可 得 到 电路 的 格林 函数 矩阵 9(w) 的 第 n 列 和 矩阵 元 ; 每 个 节点 均 操 作 一 次 , 便 可 得 到 整个 电路 的 格 书 
函数 矩阵 ， 然 后 根据 关系 Tw) = 9-1(w) 得 出 拉 普 拉 辛 矩阵 以 及 相应 的 矩阵 A 晤 引得 到 实验 测量 的 拉 普 拉 辛 
矩阵 了 (w), 就 可 以 计算 其 导 纳 谱 及 相应 的 左 / 右 本 征 矢量 , 以 及 由 此 定义 的 一 切 物理 量 , 从 而 和 理论 相 比较 . 

需要 注意 的 是 ,上 述 分 析 是 基于 电路 系统 在 AC 电流 驱动 下 的 稳 态 响应 , 即 式 (A). 而 事实 上 , 在 驱动 频率 接 


[0 


À 
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近 本 和 
实验 上 获得 稳 态 响应 的 方法 
实现 对 稳 态 响应 信和 号 的 


本 和 和 


减 


ASA 


公式 F(we)V =0 5 


可 以 求 出 本 征 频 


FE 频 率 we 一 般 也 为 复数 
(Imlwe] > 0). 对 于 L 


HIE MAR 
z Bbc 
频率 we ] 
各 应 的 抑制 , 这 里 通过 选取 合适 的 Ro 达到 此 


可 以 采用 无 磁 通 


0) 的 


界 条 件 下 , 非 互 易 AA 模型 


SUE, 待 暂 态 响应 消 


9 


FE 频 率 时 , 除了 会 产生 明显 的 共振 响应 外 , 还 会 激发 电路 系统 的 其 他 本 征 模 式 , 这 类 响应 通常 被 称 为 暂 态 响 
逝 后 再 使 用 锁 相 放大 器 (lock-in amplifier) 


的 本 征 能 量 会 出 现 复数 , 相应 地 , 此 电路 的 


图 ], 这 就 意味 着 此 频率 下 的 暂 态 响应 会 随 着 时 间 发 散 mwe] < 0) RA 


R 为 驱动 频率 的 系统 , 发 散 的 产生 不 利于 系统 响应 的 稳定 , 因此 必须 考虑 对 暂 


的 . 


Ex 


其 


ZN 


4 


FFER 


(本 


1, r= C,/Co; an KIRA 


的 Ro 以 及 C, 使 所 有 本 和 
统 以 特定 频率 w 驱动 时 , 系统 的 稳定 响应 将 以 w 模式 为 主 eH. IRD 
的 条 件 为 


另外 , 为 了 尽 可 能 多 
分 布 于 复 平 面 的 原点 附近 , 相 
y= an = 0 可 以 得 到 合适 的 好 
简洁 . 除 特别 说 明 外 , 以 1 
利用 SPICE 模拟 周期 
(10 wH, 0.4 uF,50,1.6 uF), 
的 取 值 (Cr, Cn, La, Rap) # 
上 方案 , 对 理论 相 图 [图 四 


在 复 平面 内 是 绕 原点 的 因 
模拟 得 到 的 A(®)/Co 计算 相位 O(@) = det A(D)/Co 随 © 的 变化 
果 显示 , 模拟 结果 与 理 


F 计 算 4 
边界 条 人 
J wo = 0.5 MHz 和 (r,y) = (4,1), 3 
公式 (1g) 和 ed) 


频率 的 虚 部 都 不 小 了 


r > max(Re[a,]) — 2 
岂 诱 导出 本 和 


动 频率 w = wo/V3 二 7, 这样 的 取 值 同 
WE w = wo/ VIFT 的 取 值 
F 下 含有 N = 21 个 节点 的 电路 [E 


文献 PI 里 证 明了 非 互 易 AA 模型 在 


件 


趋 肤 与 局 域 模式 的 苋 争 . 
在 电路 的 设计 上 , 只 需要 将 
界 条 件 . 


aif. 


数 下 的 竞争 关系 . KEK 


下 表现 出 趋 肤 态 , 而 在 周 


a lb) 4 


采 吻 合 得 很 好 . 


路 (此 时 的 A 与 w 无 关 ) 计算 本 和 


[频率 的 虚 部 , 从 而 估算 出 所 需要 的 Ro. 将 


[(2+r) — 28 — in] y, (22) 
Tayern-a- F) ; (23) 


ERE A/Co 的 第 n 个 本 征 值 , 在 周期 边界 条 
F 0, 即 min, (Im[w a]) > 0, 则 系统 的 响应 不 会 随时 间 发 散 . 当 系 
附录 的 推导 因 , 可 以 得 出 使 电路 响应 不 发 散 


件 下 一 般 为 复数 EI. 因此 , 只 要 选取 合适 


UD 


H y > max ( 


Telen ) (24) 


2 +r — Refan] 


a BR AG BEAL TAME 
应 地 , ERE A = 0)/Co 的 本 征 值 a 也 具有 同样 的 特点 ， 所 以 根据 式 Pa), 设 


普 中 间 ， 由 于 非 互 易 AA 模型 H 的 能 谱 


于 边界 与 周 


普 拉 溯 矩阵 以 及 相应 的 


期 边界 条 件 下 


同样 地 , 可 以 利用 与 上 一 节 周 期 边 
ERE A/Co, 然 
个 相对 


z ika) 


时 可 以 保证 公式 (5) 中 了 = R51E 尽 可 能 


ia), 基本 的 元 件 取 值 为 (Lo, Co, Ro, Cr) = 
区 动 频 率 选 为 w = wo/V6 ~ 0.2 MHz, 其 他 元 件 


由 模型 参数 确定 , 准 


同期 势 的 周期 参数 选 为 8 = 13/21. 利用 以 
! 的 三 个 典型 区 域 进行 SPICE 模拟 . 在 v = 二 1 的 拓扑 区 , 模拟 的 A/Co 本 征 谱 
PA), 而 在 缠绕 数 v = 0 的 局 域 区 ， 
图 区 qd)], 可 以 得 到 相应 缠绕 数 的 模拟 值 . 结 


I 变 成 了 实 轴 上 的 直线 (Ac) 右 图 ]; 用 


SI 


期 边界 条 件 下 的 相 图 一 致 , 只 是 拓扑 相 区 的 态 在 开 边 界 条 
为 扩展 态 . 本 节 将 同样 利用 电路 的 拉 普 拉 辛 方法 模拟 开 边 界 条 件 


V. 开 边 界 条 件 下 赵 肤 与 局 域 模式 竞争 的 模拟 


! 两 端的 开关 同时 接 入 端口 o 并 且 令 X 相关 模块 均 开 路 , 即 可 实现 开 边 
条 件 类 似 的 方法 , 通过 SPICE 模拟 重 构 出 开 边 界 条 件 下 的 电路 拉 
后 用 此 矩阵 计算 出 相应 的 本 征 左 / 右 矢 , 即 可 看 到 趋 肤 和 局 域 模式 在 不 同 参 


ITIK, 无 须 将 AC 电流 源 依次 接 入 每 个 节点 也 可 以 达到 此 目的 . 


-00008v1 


202202 


a 
a 


chinaXiv 


10 


1 n 21 0 IPR 1 


图 3. (a, b) 由 SPICE 模拟 得 到 的 电压 响应 (已 经 归 一 化 ) 在 节点 上 的 分 布 , 分 别 对 应 图 bla) 1a = 0.5 和 入 = 0.5 的 两 条 虚线 . 
频率 为 w = wo/V6 的 AC 电流 源 接 在 第 11 个 节点 上 . (c, d) 分 别 由 (a,b) 中 的 电压 分 布 根据 公式 人 计算 的 IPR, HPR 
模拟 值 , 虚线 为 理论 值 . 箭头 指 的 是 最 小 模拟 值 , 虚线 标 出 的 是 理论 相 变 值 . 

Fig. 3. (a, b) Node distributions of voltages (normalized) simulated by SPICE along dashed lines of a = 0.5 and 和 =0.5in 
Fig. 2(a), respectively. The AC current source with w = wo/V6 is connected to the 11th node. (c, d) IPRs of the voltage 
distributions in (a, b), respectively, calculated by Eq. ba), where diamonds (dashed lines) are the simulation (theoretical) 
results. The arrows indicate the minima of simulated IPRs, while the dashed lines indicate the phase transition points in 
theory. 


电路 的 格林 函数 形式 种 可 以 用 拉 普 拉 辛 的 本 征 左 / 右 矢 表示 为 
vger- E PROF o), (25) 


这 里 , a Sarat AeA AMIE AL / A RIE ECAA eS) ae BY A E 
:一 节点 接 入 AC 电流 源 , 测量 相应 的 电压 响应 , 根据 上 式 可 知 , 此 电压 响应 是 同一 相 区 下 所 有 本 征 右 矢 YE (w) 

按 系数 jr (OVE (OTİT 线性 对 加 的 结果 , 所 以 必然 表现 为 趋 肤 或 局 域 效应 . 

取 与 周期 边界 条 件 时 相同 的 元 件 参数 , 除了 7 = 0 ( 即 Ro 开路 )， 这 是 因为 开 边 界 条 件 下 所 有 的 本 征 频率 


chinaXiv:202202.00008v1 


都 是 实数 , te b,c) 所 示 , 不 存在 响应 发 散 的 问题 , 所 以 不 需要 用 电阻 抑制 发 散 . 另外 , 需要 拉 普 拉 辛 矩阵 中 的 
D(w) = 0, 否则 当 此 项 很 大 时 , 所 有 本 征 值 j(w) 趋向 于 常数 j, 由 公式 E) TA V > FOZ, 其 正比 于 输入 电 
流 , 无 法 反映 出 竞争 关系 , 这 也 是 取 y= 0 的 一 个 原因 |. 
用 SPICE 同样 模拟 了 开 边 界 条 件 下 N = 21 个 节点 的 电路 系统 , 在 节点 m = 11 接 入 频率 为 w = wo/V6 的 
AC 电流 源 , 然后 测量 每 个 节点 频率 为 w 的 电压 幅 值 , 得 到 如 图 Bla) 和 bl) 所 示 结 果 , 清楚 地 显示 在 拓扑 相 区 
时 , 响应 电压 分 布 在 右 /左边 界 附近 , 表现 为 右 / 左 趋 肤 态 ; 在 非 拓扑 相 区 时 , 响应 电压 始终 分 布 在 驱动 节点 附近 ， 
表现 为 局 域 态 . 为 了 表征 相应 电压 的 局 域 化 程度 , 定义 倒 参与 率 (inverse participation ratio, 简称 IPR): 


Pa (> val) /Ol MP): (26) 


rl Bc) #1 pla) 显示 IPR 最 低 的 点 与 理论 相 变 点 很 接近 , 因为 趋 肤 和 局 域 态 都 有 较 大 的 局 域 性 , 对 应 较 大 的 IPR 
A, 而 相 变 点 附近 扩展 性 最 强 , 对 应 的 IPR 值 都 很 小 . 在 左右 趋 肤 的 相 边界 ( 即 v = 41 之 间 的 边界 ), 对 应 于 互 
易 模型 , 因为 趋 肤 效应 消失 , 所 以 其 本 征 态 为 扩展 态 ; 在 趋 肤 与 局 域 的 相 边界 ( 即 vy = 士 1 与 v = 0 之 间 的 边界 )， 
对 应 于 趋 肤 与 局 域 竞争 的 平衡 , 同样 为 扩展 态 . 


VI. 总结 


本 文通 过 构建 经 典 电路 , 将 其 拉 普 拉 辛 矩 阵 与 非 互 易 AA 模型 的 哈密 顿 量 矩 阵 对 应 , 利用 SPICE 成 功 模拟 

THES AA 模型 重要 的 稳 态 性 质 , 包括 周期 边界 条 件 下 体现 系统 非 厄 米 拓扑 性 质 的 复 能 谱 和 能 谱 缠 绕 数 , 以 及 
开 边界 下 非 厄 米 趋 肤 效 应 与 准 无 序 局 域 化 的 竞争 ， 其 中 , 详细 讨论 了 电路 参数 的 设置 原则 和 理论 依据 , 为 进一步 
的 实验 实现 提供 了 具体 的 指导 方案 . 由 于 方案 的 普 适 性 , 本 文中 所 讨论 的 设计 原则 和 理论 可 以 直接 应 用 于 其 他 量 
子 紧 束缚 模型 的 模拟 和 实验 , 比如 文献 PA 中 提 到 的 非 互 易 AA 模型 的 对 偶 模 型 , 仅 需 将 图 四 中 的 电路 方案 适当 
修改 , 去 除 节点 间 的 INIC 元 件 以 实现 互 易 跃迁 , 并 调节 相应 的 接地 元 件 以 实现 准 周期 复 在 位 势 即 可 . 
本 文中 的 SPICE 模拟 使 用 了 LTspice 软件 . 为 了 更 接近 于 理论 结果 , 这 里 采用 的 电容 、 阻 均 为 理 
想 线性 器 件 , 并 且 将 INIC 中 放大 器 的 开 环 增益 倍数 与 输入 阻抗 分 别 设置 为 500 G 和 500 GQ 用 于 模拟 理想 放大 
器 . 由 于 实际 器 件 的 非 理想 性 , 模拟 或 者 实验 结果 可 能 会 有 一 定 偏差 , 需要 根据 情况 具体 分 析 . 
本 文 只 涉及 对 量子 紧 束 缚 模型 稳 态 性 质 的 模拟 , 实际 上 , 经 典 电路 也 可 以 用 于 对 动力 学 性 质 的 模拟 . 对 于 非 
互 易 AA 模型 , 文献 PT) 已 有 讨论 . 另外 , 由 于 电路 元 件 的 丰富 特性 , 同样 可 以 利用 非 线性 的 电路 元 件 实现 对 非 
线性 量子 系统 的 模拟 [64]. FLL, 对 于 模拟 量子 系统 而 言 , 经 典 电路 是 一 个 成 本 低廉 、 技 术 成 熟 、 模 拟 范围 广 的 
有 力 平台 . 
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附录 A: 电路 响应 不 发 散 的 条 件 


根据 本 征 频率 的 表达 式 PS), 如 使 电路 响应 不 发 散 , 需要 满足 mw /wo] > 0, 即 


0o>2m| H VQ+r) -an | = g= R OFr) tan Ne (2 十 门 十 on 


要 -1< (2-@4r +004 |V- eta] | (A1) 


这 里 将 根 式 V' 整体 看 作 复 数 , 由 于 不 等 式 左边 为 非 正 数 (7 > 0), 所 以 只 需 考虑 不 等 式 右 边 (为 根 式 的 两 倍 实 部 ) 


为 非 正 数 时 满足 的 条 件 即 可 . 因 


k 将 不 等 式 (Al) 两边 平 方 并 化 简 , 得 


一 


* 


2 * 2 2 
X +241) ee (Q+r) +a, TEA , (A2) 
两 边 再 平方 并 化 简 , 得 
(2 +r — Rela,])y* > (Im[a,])?. (A3) 
于 不 等 式 右 方 为 非 负 数 , 所 以 必须 要 求 

2+r—Rela,]>0 = r> Reļan]-—2, (A4) 
y > |Im[an]|/ V2 +r — Refan]. (A5) 

因为 需 保证 所 有 本 征 频率 都 满足 以 上 条 件 , 所 以 要 求 

> |Im[a,,]| ) 

r> max(Reļan]) 2 7 > max ( E (A6) 
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